Characterisation of laminar and vascular spatiotemporal dynamics of CBV and BOLD signals using VASO and ME-GRE at 7T in humans.

Interpretation of cortical laminar functional magnetic resonance imaging (fMRI) activity requires detailed knowledge of the spatiotemporal haemodynamic response across vascular compartments due to the well-known vascular biases (e.g. the draining veins). Further complications arise from the spatiotemporal hemodynamic response that differs depending on the duration of stimulation. This information is crucial for future studies using depth-dependent cerebral blood volume (CBV) measurements, which promise higher specificity for the cortical microvasculature than the blood oxygenation level dependent (BOLD) contrast. To date, direct information about CBV dynamics with respect to stimulus duration, cortical depth and vasculature is missing in humans. Therefore, we characterized the cortical depth-dependent CBV-haemodynamic responses across a wide set of stimulus durations with 0.9 mm isotropic spatial and 0.785 seconds effective temporal resolution in humans using slice-selective slab-inversion vascular space occupancy (SS-SI VASO). Additionally, we investigated signal contributions from macrovascular compartments using fine-scale vascular information from multi-echo gradient-echo (ME-GRE) data at 0.35 mm isotropic resolution. In total, this resulted in >7.5h of scanning per participant (n=5). We have three major findings: (I) While we could demonstrate that 1 second stimulation is viable using VASO, more than 12 seconds stimulation provides better CBV responses in terms of specificity to microvasculature, but durations beyond 24 seconds of stimulation may be wasteful for certain applications. (II) We observe that CBV responses show dilation patterns across the cortex. (III) While we found increasingly strong BOLD signal responses in vessel-dominated voxels with longer stimulation durations, we found increasingly strong CBV signal responses in vessel-dominated voxels only until 4 second stimulation durations. After 4 seconds, only the signal from non-vessel dominated voxels kept increasing. This might explain why CBV responses are more specific to the underlying neuronal activity for long stimulus durations.

Longitudinal motor system changes from acute to chronic spinal cord injury.

BACKGROUND AND PURPOSE: In acute spinal cord injury (SCI), magnetic resonance imaging (MRI) reveals tissue bridges and neurodegeneration for 2 years. This 5-year study aims to track initial lesion changes, subsequent neurodegeneration, and their impact on recovery. METHODS: This prospective longitudinal study enrolled acute SCI patients and healthy controls who were assessed clinically-and by MRI-regularly from 3 days postinjury up to 60 months. We employed histologically cross-validated quantitative MRI sequences sensitive to volume, myelin, and iron changes, thereby reflecting indirectly processes of neurodegeneration and neuroinflammation. General linear models tracked lesion and remote changes in volume, myelin- and iron-sensitive magnetic resonance indices over 5 years. Associations between lesion, degeneration, and recovery (using the Spinal Cord Independence Measure [SCIM] questionnaire and the International Standards for Neurological Classification of Spinal Cord Injury total motor score) were assessed. RESULTS: Patients' motor scores improved by an average of 12.86 (95% confidence interval [CI] = 6.70-19.00) points, and SCIM by 26.08 (95% CI = 17.00-35.20) points. Within 3-28 days post-SCI, lesion size decreased by more than two-thirds (3 days: 302.52 ± 185.80 mm2 , 28 days: 76.77 ± 88.62 mm2 ), revealing tissue bridges. Cervical cord and corticospinal tract volumes transiently increased in SCI patients by 5% and 3%, respectively, accompanied by cervical myelin decreases and iron increases. Over time, progressive atrophy was observed in both regions, which was linked to early lesion dynamics. Tissue bridges, reduced swelling, and myelin content decreases were predictive of long-term motor score recovery and improved SCIM score. CONCLUSIONS: Studying acute changes and their impact on longer follow-up provides insights into SCI trajectory, highlighting the importance of acute intervention while indicating the potential to influence outcomes in the later stages.

Fiber-orientation independent component of R2* obtained from single-orientation MRI measurements in simulations and a post-mortem human optic chiasm.

The effective transverse relaxation rate (R2*) is sensitive to the microstructure of the human brain like the g-ratio which characterises the relative myelination of axons. However, the fibre-orientation dependence of R2* degrades its reproducibility and any microstructural derivative measure. To estimate its orientation-independent part (R2,iso*) from single multi-echo gradient-recalled-echo (meGRE) measurements at arbitrary orientations, a second-order polynomial in time model (hereafter M2) can be used. Its linear time-dependent parameter, ß1, can be biophysically related to R2,iso* when neglecting the myelin water (MW) signal in the hollow cylinder fibre model (HCFM). Here, we examined the performance of M2 using experimental and simulated data with variable g-ratio and fibre dispersion. We found that the fitted ß1 can estimate R2,iso* using meGRE with long maximum-echo time (TEmax ≈ 54 ms), but not accurately captures its microscopic dependence on the g-ratio (error 84%). We proposed a new heuristic expression for ß1 that reduced the error to 12% for ex vivo compartmental R2 values. Using the new expression, we could estimate an MW fraction of 0.14 for fibres with negligible dispersion in a fixed human optic chiasm for the ex vivo compartmental R2 values but not for the in vivo values. M2 and the HCFM-based simulations failed to explain the measured R2*-orientation-dependence around the magic angle for a typical in vivo meGRE protocol (with TEmax ≈ 18 ms). In conclusion, further validation and the development of movement-robust in vivo meGRE protocols with TEmax ≈ 54 ms are required before M2 can be used to estimate R2,iso* in subjects.

Dissociating distinct cortical networks associated with subregions of the human medial temporal lobe using precision neuroimaging.

Tract-tracing studies in primates indicate that different subregions of the medial temporal lobe (MTL) are connected with multiple brain regions. However, no clear framework defining the distributed anatomy associated with the human MTL exists. This gap in knowledge originates in notoriously low MRI data quality in the anterior human MTL and in group-level blurring of idiosyncratic anatomy between adjacent brain regions, such as entorhinal and perirhinal cortices, and parahippocampal areas TH/TF. Using MRI, we intensively scanned four human individuals and collected whole-brain data with unprecedented MTL signal quality. Following detailed exploration of cortical networks associated with MTL subregions within each individual, we discovered three biologically meaningful networks associated with the entorhinal cortex, perirhinal cortex, and parahippocampal area TH, respectively. Our findings define the anatomical constraints within which human mnemonic functions must operate and are insightful for examining the evolutionary trajectory of the MTL connectivity across species.

Reliability of resting-state functional connectivity in the human spinal cord: Assessing the impact of distinct noise sources.

The investigation of spontaneous fluctuations of the blood-oxygen-level-dependent (BOLD) signal has recently been extended from the brain to the spinal cord, where it has stimulated interest from a clinical perspective. A number of resting-state functional magnetic resonance imaging (fMRI) studies have demonstrated robust functional connectivity between the time series of BOLD fluctuations in bilateral dorsal horns and between those in bilateral ventral horns, in line with the functional neuroanatomy of the spinal cord. A necessary step prior to extension to clinical studies is assessing the reliability of such resting-state signals, which we aimed to do here in a group of 45 healthy young adults at the clinically prevalent field strength of 3T. When investigating connectivity in the entire cervical spinal cord, we observed fair to good reliability for dorsal-dorsal and ventral-ventral connectivity, whereas reliability was poor for within- and between-hemicord dorsal-ventral connectivity. Considering how prone spinal cord fMRI is to noise, we extensively investigated the impact of distinct noise sources and made two crucial observations: removal of physiological noise led to a reduction in functional connectivity strength and reliability - due to the removal of stable and participant-specific noise patterns - whereas removal of thermal noise considerably increased the detectability of functional connectivity without a clear influence on reliability. Finally, we also assessed connectivity within spinal cord segments and observed that while the pattern of connectivity was similar to that of whole cervical cord, reliability at the level of single segments was consistently poor. Taken together, our results demonstrate the presence of reliable resting-state functional connectivity in the human spinal cord even after thoroughly accounting for physiological and thermal noise, but at the same time urge caution if focal changes in connectivity (e.g. due to segmental lesions) are to be studied, especially in a longitudinal manner.

Coherent, time-shifted patterns of microstructural plasticity during motor-skill learning.

Motor skill learning relies on neural plasticity in the motor and limbic systems. However, the spatial and temporal characteristics of these changes-and their microstructural underpinnings-remain unclear. Eighteen healthy males received 1 h of training in a computer-based motion game, 4 times a week, for 4 consecutive weeks, while 14 untrained participants underwent scanning only. Performance improvements were observed in all trained participants. Serial myelin- and iron-sensitive multiparametric mapping at 3T during this period of intensive motor skill acquisition revealed temporally and spatially distributed, performance-related microstructural changes in the grey and white matter across a corticospinal-cerebellar-hippocampal circuit. Analysis of the trajectory of these transient changes suggested time-shifted cascades of plasticity from the dominant sensorimotor system to the contralateral hippocampus. In the cranial corticospinal tracts, changes in myelin-sensitive metrics during training in the posterior limb of the internal capsule were of greater magnitude in those who trained their upper limbs vs. lower limb trainees. Motor skill learning is associated with waves of grey and white matter plasticity, across a broad sensorimotor network.

Germany's journey toward 14 Tesla human magnetic resonance.

Multiple sites within Germany operate human MRI systems with magnetic fields either at 7 Tesla or 9.4 Tesla. In 2013, these sites formed a network to facilitate and harmonize the research being conducted at the different sites and make this technology available to a larger community of researchers and clinicians not only within Germany, but also worldwide. The German Ultrahigh Field Imaging (GUFI) network has defined a strategic goal to establish a 14 Tesla whole-body human MRI system as a national research resource in Germany as the next progression in magnetic field strength. This paper summarizes the history of this initiative, the current status, the motivation for pursuing MR imaging and spectroscopy at such a high magnetic field strength, and the technical and funding challenges involved. It focuses on the scientific and science policy process from the perspective in Germany, and is not intended to be a comprehensive systematic review of the benefits and technical challenges of higher field strengths.

Reliability of spinal cord measures based on synthetic T1-weighted MRI derived from multiparametric mapping (MPM).

Short MRI acquisition time, high signal-to-noise ratio, and high reliability are crucial for image quality when scanning healthy volunteers and patients. Cross-sectional cervical cord area (CSA) has been suggested as a marker of neurodegeneration and potential outcome measure in clinical trials and is conventionally measured on T1-weigthed 3D Magnetization Prepared Rapid Acquisition Gradient-Echo (MPRAGE) images. This study aims to reduce the acquisition time for the comprehensive assessment of the spinal cord, which is typically based on MPRAGE for morphometry and multi-parameter mapping (MPM) for microstructure. The MPRAGE is replaced by a synthetic T1-w MRI (synT1-w) estimated from the MPM, in order to measure CSA. SynT1-w images were reconstructed using the MPRAGE signal equation based on quantitative maps of proton density (PD), longitudinal (R1) and effective transverse (R2*) relaxation rates. The reliability of CSA measurements from synT1-w images was determined within a multi-center test-retest study format and validated against acquired MPRAGE scans by assessing the agreement between both methods. The response to pathological changes was tested by longitudinally measuring spinal cord atrophy following spinal cord injury (SCI) for synT1-w and MPRAGE using linear mixed effect models. CSA measurements based on the synT1-w MRI showed high intra-site (Coefficient of variation [CoV]: 1.43% to 2.71%) and inter-site repeatability (CoV: 2.90% to 5.76%), and only a minor deviation of -1.65 mm2 compared to MPRAGE. Crucially, by assessing atrophy rates and by comparing SCI patients with healthy controls longitudinally, differences between synT1-w and MPRAGE were negligible. These results demonstrate that reliable estimates of CSA can be obtained from synT1-w images, thereby reducing scan time significantly.

High-resolution quantitative and functional MRI indicate lower myelination of thin and thick stripes in human secondary visual cortex.

The characterization of cortical myelination is essential for the study of structure-function relationships in the human brain. However, knowledge about cortical myelination is largely based on post-mortem histology, which generally renders direct comparison to function impossible. The repeating pattern of pale-thin-pale-thick stripes of cytochrome oxidase (CO) activity in the primate secondary visual cortex (V2) is a prominent columnar system, in which histology also indicates different myelination of thin/thick versus pale stripes. We used quantitative magnetic resonance imaging (qMRI) in conjunction with functional magnetic resonance imaging (fMRI) at ultra-high field strength (7 T) to localize and study myelination of stripes in four human participants at sub-millimeter resolution in vivo. Thin and thick stripes were functionally localized by exploiting their sensitivity to color and binocular disparity, respectively. Resulting functional activation maps showed robust stripe patterns in V2 which enabled further comparison of quantitative relaxation parameters between stripe types. Thereby, we found lower longitudinal relaxation rates (R1) of thin and thick stripes compared to surrounding gray matter in the order of 1-2%, indicating higher myelination of pale stripes. No consistent differences were found for effective transverse relaxation rates (R2*). The study demonstrates the feasibility to investigate structure-function relationships in living humans within one cortical area at the level of columnar systems using qMRI.

Dynamics of progressive degeneration of major spinal pathways following spinal cord injury: A longitudinal study.

BACKGROUND: Following spinal cord injury (SCI), disease processes spread gradually along the spinal cord forming a spatial gradient with most pronounced changes located at the lesion site. However, the dynamics of this gradient in SCI patients is not established. OBJECTIVE: This study tracks the spatiotemporal dynamics of remote anterograde and retrograde spinal tract degeneration in the upper cervical cord following SCI over two years utilizing quantitative MRI. METHODS: Twenty-three acute SCI patients (11 paraplegics, 12 tetraplegics) and 21 healthy controls were scanned with a T1-weighted sequence for volumetry and a FLASH sequence for myelin-sensitive magnetization transfer saturation (MTsat) of the upper cervical cord. We estimated myelin content from MTsat maps within the corticospinal tracts (CST) and dorsal columns (DC) and measured spinal cord atrophy by means of left-right width (LRW) and anterior-posterior width (APW) on the T1-weighted images across cervical levels C1-C3. MTsat in the CST and LRW were considered proxies for retrograde degeneration, while MTsat in the DC and APW provided evidence for anterograde degeneration, respectively. Using regression models, we compared the temporal and spatial trajectories of these MRI readouts between tetraplegics, paraplegics, and controls over a 2-year period and assessed their associations with clinical improvement. RESULTS: Linear rates and absolute differences in myelin-sensitive MTsat indicated retrograde and anterograde neurodegeneration in the CST and DC, respectively. Changes in MTsat within the CST and in LRW progressively developed over time forming a gradient towards lower cervical levels by 2 years after injury, especially in tetraplegics (change per cervical level in MTsat: -0.247 p.u./level, p = 0.034; in LRW: -0.323 mm/level, p = 0.024). MTsat within the DC was already decreased at cervical levels C1-C3 at baseline (1.5 months after injury) in both tetra- and paraplegics, while linear decreases in APW over time were similar across C1-C3, preserving the spatial gradient. The relative improvement in light touch score was associated with MTsat within the DC at baseline (rs = 0.575, p = 0.014). CONCLUSION: Rostral and remote to the injury, the CST and DC show ongoing structural changes, indicative of myelin reductions and atrophy within 2 years after SCI. While anterograde degeneration in the DC was already detectable uniformly at C1-C3 early following SCI, retrograde degeneration in the CST developed over time revealing specific spatial and temporal neurodegenerative gradients. Disentangling and quantifying such dynamic pathological processes may provide biomarkers for regenerative and remyelinating therapies along entire spinal pathways.

HiHi fMRI: a data-reordering method for measuring the hemodynamic response of the brain with high temporal resolution and high SNR.

There is emerging evidence that sampling the blood-oxygen-level-dependent (BOLD) response with high temporal resolution opens up new avenues to study the in vivo functioning of the human brain with functional magnetic resonance imaging. Because the speed of sampling and the signal level are intrinsically connected in magnetic resonance imaging via the T1 relaxation time, optimization efforts usually must make a trade-off to increase the temporal sampling rate at the cost of the signal level. We present a method, which combines a sparse event-related stimulus paradigm with subsequent data reshuffling to achieve high temporal resolution while maintaining high signal levels (HiHi). The proof-of-principle is presented by separately measuring the single-voxel time course of the BOLD response in both the primary visual and primary motor cortices with 100-ms temporal resolution.

Quantitative MRI maps of human neocortex explored using cell type-specific gene expression analysis.

Quantitative magnetic resonance imaging (qMRI) allows extraction of reproducible and robust parameter maps. However, the connection to underlying biological substrates remains murky, especially in the complex, densely packed cortex. We investigated associations in human neocortex between qMRI parameters and neocortical cell types by comparing the spatial distribution of the qMRI parameters longitudinal relaxation rate (${R_{1}}$), effective transverse relaxation rate (${R_{2}}^{\ast }$), and magnetization transfer saturation (MTsat) to gene expression from the Allen Human Brain Atlas, then combining this with lists of genes enriched in specific cell types found in the human brain. As qMRI parameters are magnetic field strength-dependent, the analysis was performed on MRI data at 3T and 7T. All qMRI parameters significantly covaried with genes enriched in GABA- and glutamatergic neurons, i.e. they were associated with cytoarchitecture. The qMRI parameters also significantly covaried with the distribution of genes enriched in astrocytes (${R_{2}}^{\ast }$ at 3T, ${R_{1}}$ at 7T), endothelial cells (${R_{1}}$ and MTsat at 3T), microglia (${R_{1}}$ and MTsat at 3T, ${R_{1}}$ at 7T), and oligodendrocytes and oligodendrocyte precursor cells (${R_{1}}$ at 7T). These results advance the potential use of qMRI parameters as biomarkers for specific cell types.

B 1 + $$ {B}_1

PURPOSE: Magnetization transfer saturation ( MTsat $$ \mathrm{MTsat} $$ ) is a useful marker to probe tissue macromolecular content and myelination in the brain. The increased B 1 + $$ {B}_1^{+} $$ -inhomogeneity at ≥ 7 $$ \ge 7 $$ T and significantly larger saturation pulse flip angles which are often used for postmortem studies exceed the limits where previous MTsat $$ \mathrm{MTsat} $$ B 1 + $$ {B}_1^{+} $$ correction methods are applicable. Here, we develop a calibration-based correction model and procedure, and validate and evaluate it in postmortem 7T data of whole chimpanzee brains. THEORY: The B 1 + $$ {B}_1^{+} $$ dependence of MTsat $$ \mathrm{MTsat} $$ was investigated by varying the off-resonance saturation pulse flip angle. For the range of saturation pulse flip angles applied in typical experiments on postmortem tissue, the dependence was close to linear. A linear model with a single calibration constant C $$ C $$ is proposed to correct bias in MTsat $$ \mathrm{MTsat} $$ by mapping it to the reference value of the saturation pulse flip angle. METHODS: C $$ C $$ was estimated voxel-wise in five postmortem chimpanzee brains. "Individual-based global parameters" were obtained by calculating the mean C $$ C $$ within individual specimen brains and "group-based global parameters" by calculating the means of the individual-based global parameters across the five brains. RESULTS: The linear calibration model described the data well, though C $$ C $$ was not entirely independent of the underlying tissue and B 1 + $$ {B}_1^{+} $$ . Individual-based correction parameters and a group-based global correction parameter ( C = 1 . 2 $$ C=1.2 $$ ) led to visible, quantifiable reductions of B 1 + $$ {B}_1^{+} $$ -biases in high-resolution MTsat $$ \mathrm{MTsat} $$ maps. CONCLUSION: The presented model and calibration approach effectively corrects for B 1 + $$ {B}_1^{+} $$ inhomogeneities in postmortem 7T data.

Brain structure and function: a multidisciplinary pipeline to study hominoid brain evolution.

To decipher the evolution of the hominoid brain and its functions, it is essential to conduct comparative studies in primates, including our closest living relatives. However, strong ethical concerns preclude in vivo neuroimaging of great apes. We propose a responsible and multidisciplinary alternative approach that links behavior to brain anatomy in non-human primates from diverse ecological backgrounds. The brains of primates observed in the wild or in captivity are extracted and fixed shortly after natural death, and then studied using advanced MRI neuroimaging and histology to reveal macro- and microstructures. By linking detailed neuroanatomy with observed behavior within and across primate species, our approach provides new perspectives on brain evolution. Combined with endocranial brain imprints extracted from computed tomographic scans of the skulls these data provide a framework for decoding evolutionary changes in hominin fossils. This approach is poised to become a key resource for investigating the evolution and functional differentiation of hominoid brains.

Reliability of task-based fMRI in the dorsal horn of the human spinal cord.

The application of functional magnetic resonance imaging (fMRI) to the human spinal cord is still a relatively small field of research and faces many challenges. Here we aimed to probe the limitations of task-based spinal fMRI at 3T by investigating the reliability of spinal cord blood oxygen level dependent (BOLD) responses to repeated nociceptive stimulation across two consecutive days in 40 healthy volunteers. We assessed the test-retest reliability of subjective ratings, autonomic responses, and spinal cord BOLD responses to short heat pain stimuli (1s duration) using the intraclass correlation coefficient (ICC). At the group level, we observed robust autonomic responses as well as spatially specific spinal cord BOLD responses at the expected location, but no spatial overlap in BOLD response patterns across days. While autonomic indicators of pain processing showed good-to-excellent reliability, both ß-estimates and z-scores of task-related BOLD responses showed poor reliability across days in the target region (gray matter of the ipsilateral dorsal horn). When taking into account the sensitivity of gradient-echo echo planar imaging (GE-EPI) to draining vein signals by including the venous plexus in the analysis, we observed BOLD responses with good reliability across days. Taken together, these results demonstrate that heat pain stimuli as short as one second are able to evoke a robust and spatially specific BOLD response, which is however strongly variable within participants across time, resulting in low reliability in the dorsal horn gray matter. Further improvements in data acquisition and analysis techniques are thus necessary before event-related spinal cord fMRI as used here can be reliably employed in longitudinal designs or clinical settings.

Reliability of multi-parameter mapping (MPM) in the cervical cord: A multi-center multi-vendor quantitative MRI study.

MRI based multicenter studies which target neurological pathologies affecting the spinal cord and brain - including spinal cord injury (SCI) - require standardized acquisition protocols and image processing methods. We have optimized and applied a multi-parameter mapping (MPM) protocol that simultaneously covers the brain and the cervical cord within a traveling heads study across six clinical centers (Leutritz et al., 2020). The MPM protocol includes quantitative maps (magnetization transfer saturation (MT), proton density (PD), longitudinal (R1), and effective transverse (R2*) relaxation rates) sensitive to myelination, water content, iron concentration, and morphometric measures, such as cross-sectional cord area. Previously, we assessed the repeatability and reproducibility of the brain MPM data acquired in the five healthy participants who underwent two scan-rescans (Leutritz et al., 2020). This study focuses on the cervical cord MPM data derived from the same acquisitions to determine its repeatability and reproducibility in the cervical cord. MPM matrices of the cervical cord were generated and processed using the hMRI and the spinal cord toolbox. To determine reliability of the cervical MPM data, the intra-site (i.e., scan-rescan) coefficient of variation (CoV), inter-site CoV, and bias within region of interests (C1, C2 and C3 levels) were determined. The range of the mean intra- and inter-site CoV of MT, R1 and PD was between 2.5% and 12%, and between 1.1% and 4.0% for the morphometric measures. In conclusion, the cervical MPM data showed a high repeatability and reproducibility for key imaging biomarkers and hence can be employed as a standardized tool in multi-center studies, including clinical trials.

Automated slice-specific z-shimming for functional magnetic resonance imaging of the human spinal cord.

Functional magnetic resonance imaging (fMRI) of the human spinal cord faces many challenges, such as signal loss due to local magnetic field inhomogeneities. This issue can be addressed with slice-specific z-shimming, which compensates for the dephasing effect of the inhomogeneities using a slice-specific gradient pulse. Here, we aim to address outstanding issues regarding this technique by evaluating its effects on several aspects that are directly relevant for spinal fMRI and by developing two automated procedures in order to improve upon the time-consuming and subjective nature of manual selection of z-shims: one procedure finds the z-shim that maximizes signal intensity in each slice of an EPI reference-scan and the other finds the through-slice field inhomogeneity for each EPI-slice in field map data and calculates the required compensation gradient moment. We demonstrate that the beneficial effects of z-shimming are apparent across different echo times, hold true for both the dorsal and ventral horn, and are also apparent in the temporal signal-to-noise ratio (tSNR) of EPI time-series data. Both of our automated approaches were faster than the manual approach, lead to significant improvements in gray matter tSNR compared to no z-shimming and resulted in beneficial effects that were stable across time. While the field-map-based approach performed slightly worse than the manual approach, the EPI-based approach performed as well as the manual one and was furthermore validated on an external corticospinal data-set (N > 100). Together, automated z-shimming may improve the data quality of future spinal fMRI studies and lead to increased reproducibility in longitudinal studies.

Error quantification in multi-parameter mapping facilitates robust estimation and enhanced group level sensitivity.

Multi-Parameter Mapping (MPM) is a comprehensive quantitative neuroimaging protocol that enables estimation of four physical parameters (longitudinal and effective transverse relaxation rates R1 and R2*, proton density PD, and magnetization transfer saturation MTsat) that are sensitive to microstructural tissue properties such as iron and myelin content. Their capability to reveal microstructural brain differences, however, is tightly bound to controlling random noise and artefacts (e.g. caused by head motion) in the signal. Here, we introduced a method to estimate the local error of PD, R1, and MTsat maps that captures both noise and artefacts on a routine basis without requiring additional data. To investigate the method's sensitivity to random noise, we calculated the model-based signal-to-noise ratio (mSNR) and showed in measurements and simulations that it correlated linearly with an experimental raw-image-based SNR map. We found that the mSNR varied with MPM protocols, magnetic field strength (3T vs. 7T) and MPM parameters: it halved from PD to R1 and decreased from PD to MTsat by a factor of 3-4. Exploring the artefact-sensitivity of the error maps, we generated robust MPM parameters using two successive acquisitions of each contrast and the acquisition-specific errors to down-weight erroneous regions. The resulting robust MPM parameters showed reduced variability at the group level as compared to their single-repeat or averaged counterparts. The error and mSNR maps may better inform power-calculations by accounting for local data quality variations across measurements. Code to compute the mSNR maps and robustly combined MPM maps is available in the open-source hMRI toolbox.

A unified 3D map of microscopic architecture and MRI of the human brain.

We present the first three-dimensional (3D) concordance maps of cyto- and fiber architecture of the human brain, combining histology, immunohistochemistry, and 7-T quantitative magnetic resonance imaging (MRI), in two individual specimens. These 3D maps each integrate data from approximately 800 microscopy sections per brain, showing neuronal and glial cell bodies, nerve fibers, and interneuronal populations, as well as ultrahigh-field quantitative MRI, all coaligned at the 200-µm scale to the stacked blockface images obtained during sectioning. These unprecedented 3D multimodal datasets are shared without any restrictions and provide a unique resource for the joint study of cell and fiber architecture of the brain, detailed anatomical atlasing, or modeling of the microscopic underpinnings of MRI contrasts.